Dynamical excitation of space-time modes of compact objects 
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We discuss, in the perturbative regime, the scattering of Gaussian pulses of odd-parity gravita- 
tional radiation off a non-rotating relativistic star and a Schwarzschild Black Hole. We focus on the 
excitation of the w-modes of the star as a function of the width b of the pulse and we contrast it 
with the outcome of a Schwarzschild Black Hole of the same mass. For sufficiently narrow values 
of 6, the waveforms are dominated by characteristic space-time modes. On the other hand, for 
sufficiently large values of b the backscattered signal is dominated by the tail of the Regge- Wheeler 
potential, the quasi-normal modes are not excited and the nature of the central object cannot be 
established. We view this work as a useful contribution to the comparison between perturbative 
results and forthcoming w-mode 3D-nonlinear numerical simulation. 

PACS numbers: 04.30.Db, 04.40.Dg, 95.30.Sf, 
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I. INTRODUCTION 

The pioneering works of Vishveshwara [ij , Press fl| and 
Davis, Ruffini and Tiomno unambiguously showed 
that a non-spherical gravitational perturbation of a 
Schwarzschild Black Hole is radiated away via exponen- 
tially damped harmonic oscillations. These oscillations 
are interpreted as space-time vibrational modes. The 
properties of these quasi-normal modes (QNMs hence- 
forth) of Black Holes have been thoroughly studied since 
then (see for example Refs. [1, H, [y] and references 
therein). Relativistic stars can also have space-time vi- 
brational modes, the so-called w-modes Q. These modes 
are purely relativistic and, contrary to fluid modes, are 
absent in Newtonian theory. The fundamental w-mode 
frequency of a typical neutron star of radius ~ 10 km 
and mass 1.4Mq is expected to lie in the range of 10 12 
kHz and to have a damping time of ~ 10~* s 0. 

The issue of the excitation of w-modes in astrophys- 
ically motivated scenarios has been deeply investigated 
in the literature. Andersson and Kokkotas [9| showed 
that, in the odd-parity case, the scattering of a Gaussian 
pulse of gravitational waves off a constant density non 
rotating star generates a waveform that, in close anal- 
ogy with the Black Hole case, is characterized by three 
phases: (i) a precursor, mainly related to the choice of 
the initial data and determined by the backscattering of 
the background curvature while the pulse is entering in 
the gravitational field of the star; (ii) a burst; (iii) a ring- 
down phase dominated by w-modes, whose presence was 
inferred by looking at the Fourier spectrum of the signals. 
Since the star is non-rotating, the signal eventually dies 
out with a power-law tail typical of Schwarzschild space- 
time [M EM- Allen and coworkers [l2| and Ruoff [H 
addressed, by means of time-domain perturbative anal- 
ysis, the same problem in the even-parity case, focusing 
on gravitational wave scattering scenarios. They consid- 
ered a large sample of initial configurations as well as 



star models of different compaction. Their main findings 
were: (i) w-modes are present only for non-conformally 
flat initial data (i.e., some radiative field needs to be 
injected in the system) and (ii) the strength of the w- 
mode signal depends on the compaction of the star. 
These pioneerin g st udies were later extended or refined 
in Refs. [13, ElM p. Ill 111 [M [111, [13, [11, [13. In 
particular, Refs. [ll, [l^ [l3| considered the scatter- 
ing off the star of particles moving along open orbits and 
realized that the w-mode excitation strongly depends on 
the orbital parameters: the closer the turning point of 
the orbit is to the star (i.e., the higher is the frequency 
of the gravitational wave instantaneously emitted by the 
particle), the larger is the presence of w-modes. Consis- 
tently, Ref. [1H| showed that (modulo a simplified treat- 
ment of the star surface) if the source of perturbation 
is a spatially extended axisymmetric distribution of fluid 
matter (like a quadrupolar shell) plunging on the star, 
the w-modes are not excited, but the energy spectrum is 
dominated by low-frequency contributions due to curva- 
ture backscattering. In addition, Ref. i26|] addressed the 
late-time decay of the trapped mode for ultra-compact, 
highly relativistic constant density stars. The presence 
of trapped w-modes in stars with a first-order phase 
transition (a density discontinuity) was also discussed in 
Ref. [13- 

In this work we analyze the problem of w-modes ex- 
citation in relativistic stars (in the perturbative regime) 
by emphasizing the analogies with the Black Hole case. 
For a given odd-parity gravitational wave multipole i, 
we consider the scattering of Gaussian pulses of gravita- 
tional radiation of different width b off relativistic stars 
(either with constant density or with a polytropic equa- 
tion of state) of 1.4 Mq and we contrast such signals with 
those emitted by a non-rotating Black Hole of the same 
mass. We focus on the excitation of space-time modes as 
a function of the width b of the Gaussian. We find that 
space-times modes can be clearly identified only if the 
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Gaussian wave-packet is sufBciently narrow (i.e., small 
6) . On the other hand, for large wave-packets the internal 
structure of the object is unaffected by the perturbation, 
tail effects are dominating and the gravitational wave- 
forms generated by stars or Black Holes are practically 
identical. 



II. NUMERICAL FRAMEWORK 

A. Relativistic stars 

From the spherically symmetric line element in 
Schwarzschild coordinates 

ds^ = -e^°'dt^ + e^^dr"^ + {d6^ + sin^ Odip^) , (2.1) 

by assuming the stress energy tensor of a perfect fluid as 
T^^^ = (p -\- ijl)u^u'^ -\- pg^'^ , where p is the pressure and 
fi the total energy density of the star, the Einstein equa- 
tions reduce to the Tolman-Oppenheimer-Volkoff (TOV) 
equations of stellar equilibrium: 



da ^ (m + Airr^p) 
dr (r^ — 2mr) ' 
dp , .da 

/--(p + M):r- (2.4) 

dr dr 

Since we are using Schwarzschild coordinates, we also 
have that e~'^^ = 1 — 2m(r)/r, where m{r) is the mass 
contained in a sphere of radius r. This system of equa- 
tions needs, to be solved, the specification of an Equation 
of State (EoS). We use the simplest two: the constant 
density EoS and an adiabatic EoS in the form p = K^^ . 
We consider two polytropic models (named A and B) and 
two constant energy density models (Ac and Be) with 
the same compactness. All the models share the same 
mass AI = IAMq, and their specific properties are listed 
in Table |T1 If not differently stated, we use geometrized 
units c = G = 1 with Mq = 1 . 

B. Odd-parity perturbations 

Odd-parity linear perturbations of Black Holes and 
neutron stars (in the absence of external matter source) 
are described by a simple linear equation 

a^^-M _ + = , (2.5) 

for a master function 5''^°) = ^f^"-* that is related to the 
metric degrees of freedom (see for example Ref. [11]). For 
a star of radius R, the potential is given by 



TABLE I; From up to down the rows report: the polytropic 
constant K, the adiabatic index F, the mass of the star M, 
its radius R, the central pressure pc, the central total energy 
density fic and the compaction parameter M/R, for all the 
stellar models considered. 



EoS 


A 


B 


Ac 


Be 


K 


56.16 


82.69 






V 


2 


2 






M 


1.40 


1.40 


1.40 


1.40 


R 


6.64 


9.10 


6.64 


9.10 


Pc 


8.84 X 10"'' 


1.83 X 10"* 


2.13 X lO""* 


4.98 X 10"^ 


He 


3.96 X 10"^ 


1.48 X 10"^ 


1.14 X 10"^ 


4.43 X 10"" 


M/R 


0.21 


0.15 


0.21 


0.15 



which reduces to the standard Regge- Wheeler poten- 
tial for r > R, where M = m{R) is the total mass 
and p — fi — 0. The latter holds also for the Black 
Hole case. We expressed Eq. (|2.5p using a tortoise 
coordinate defined as dr^/dr = exp(/3 — a). This re- 
duces in vacuum to the Regge- Wheeler tortoise coordi- 
nate r* = r -|- 2Mln(r/(2M) — 1). The power emitted in 
gravitational waves is: 

^(o)^^^(o)^^|±2)!|^(„)|,^ (2.7) 
e>2 e>2 ^ 

where the over-dot stands for coordinate time derivative. 



C. Initial data and simulation method 

For NS and Black Holes Eq. (|2.5p is solved in the time 
domain as an initial value problem. In the case of the 
polytropic EOS, we need first to integrate numerically 
the TOV equations (|2.2p - (|2.4[) to compute the potential 
V}°\ For a given central pressure Pc (see Table |T]) the 
TOV equations are integrated numerically (from the cen- 
ter outward) using a standard fourth-order Runge-Kutta 
integration scheme with adaptive step size. 

As initial data for (\E'^°\ 9(5'^°^) we set up an ingoing 
(5j\I/(°) = 9r.*'°^) Gaussian pulse of tunable width b 

=AAexp[-(r-ro)V&'] , (2.8) 

where A/" is a normalization constant determined by 
equating to one the integral of Eq. (|2.8p all over the ra- 
dial domain. This is a simple, but sufficiently general, 
way to represent a "distortion" of the space-time (whose 
intimate origin depends on the particular astrophysical 
setting), and to introduce in the system a proper scale 
through the width of the Gaussian. 

Let us also summarize the basilar elements of our nu- 
merical procedure. Eq. (|2.5p has been discretized on an 
evenly spaced grid (in for the Black Hole and in r for 
the star) and solved using a standard implementation of 
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the second-order Lax-WendrofT method as implemented 
for example in [25] . We have performed convergence tests 
of the code which assured a convergence factor of ~ 2. 
A resolution of Ar* — 0.01 and Ar = 0.015 is sufficient 
to be in the convergence regime. Since we have imple- 
mented standard Sommerfeld outgoing boundary condi- 
tions (see Ref. [s^ for improved, non-reflecting boundary 
conditions), we can't avoid some spurious reflections to 
come back from boundaries. To avoid that this effect 
contaminates too much the late-time tails of the signals, 
we need to choose radial grids sufficiently extended, say 
n e [-2000, 2000] and r e [0, 2000]. 



III. RESULTS 



A. Analysis of the waveforms 




150 200 250 300 350 400 



We analyzed the gravitational wave response of rela- 
tivistic stars described by the four (two polytropic and 
two constant energy density) models in Tab. |T] and of a 
Black Hole of the same mass to an impinging gravita- 
tional wave-packet of the form (|2.8p . We focus on the 
dependence of the excitation of the star w-modes (and of 
the Black Hole QNMs) ring-down on the width b. The 
Gaussian is centered at rg — 100; the waveforms are ex- 
tracted at r°^^ = 900 {r°^^ ~ 916) and shown versus 
observer retarded time u — t — r°^^. Fig. [1] exhibits the 
waveforms, for Model A, Model B and the black hole for 
b = 2 (top), & = 8 (middle) and 6 = 20 (bottom). The 
main panel depicts the modulus on a logarithmic scale, 
in order to highlight the late-time non-oscillatory tail. 

Let us first discuss the main features of the signal of 
Fig. [U starting with the "narrow" pulse, b = 2. In the 
case of the Black Hole, the ring-down has the "standard" 
shape dominated by the fundamental mode that is quoted 
in textbooks. In the case of the stars, a damped harmonic 
oscillations due to w-modes appears (we shall make this 
statement more precise below) . The waveforms show the 
common global behavior precursor - burst - ring-down - 
tail. The precursor is determined by the choice of ini- 
tial data and by the long-range features of the potential; 
this implies that, until u ~ 100, the three waveforms 
are superposed. At later times, the short-range struc- 
ture (burst-ring-down) becomes apparent. For the Black 
Hole the burst is related to the pulse passing through 
the peak of Regge- Wheeler potential. After the pulse 
the quasi-harmonic oscillatory regime shows up. When 
b is increased {b = 8), the features remain unchanged, 
but, although the non-oscillatory tail is not dominating 
yet, the amplitude of the damped oscillation is smaller 
and lasting for a shorter time. A further enlargement of 
the Gaussian causes the ingoing pulse to be almost com- 
pletely reflected back by the "tail" of the potential, so 
that the emerging waveform is unaffected by the prop- 
erties of the central object. The bottom panel of Fig. [1] 
highlights this effect for b = 20: no quasi-normal oscil- 
lations are present. It turns out that the waveforms are 
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FIG. 1: Dependence of the ring-down phase on the width 6 of 
the Gaussian pulse: for b = 2 {top panel) and & = 8 {middle 
panel) the process of excitation of the space-time modes shows 
the same qualitative features for the Black Hole and for the 
star. The waveforms for 6 = 20 {bottom panel) show there is 
basically no difference between the gravitational wave signal 
backscattered from a stars of TableUand from a Schwarzschild 
Black Hole with the same mass. 
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FIG. 3: The ratio, as a function of b, of the energy released 
in gravitational waves by the Black Hole and different star 
models. 



FIG. 2: Energy spectra (from Model A) for different values of 
b. The maximum frequency is consistent with tJmax — S\/2/b. 
See text for discussion. 



perfectly superposed and any characteristic signature of 
the Black Hole or of the star (for any star model, see 
below) disappears. We have checked through a linear 
fit that the tail is (asymptotically) in perfect agreement 
with the Price law: t-^'^+^ [13, [III- 

The absence of QNMs for large values of b is qualita- 
tively explained by means of the following argument (see 
also Sec. IX of Ref. [1^): in the frequency domain, the 
Gaussian perturbation Eq. (|2.8p is equivalent to a Gaus- 
sian of variance a^j = V^/b and contains all frequencies. 
However this means that the amplitudes of the modes 
excited by this kind of initial data will be exponentially 
suppressed if their frequencies are greater than the one 
corresponding to three standard deviations, i.e., if their 
frequency is greater than a sort of effective maximum fre- 
quency given by w^ax — ^Ci^; — SV^/b. Generally speak- 
ing, we expect to trigger the space-time modes of the 
star (or of the Black Hole) only when b is such that w^^x 
is larger than the frequency of the least damped quasi- 
normal mode of the system. In order to show how this 
argument works, let us note that we have Wj^-^^x — 2.12, 
w^ax - 0-53, t^rnix - 0-30 and uj^^^ ~ 0.21. TablelUlists 
the first six w-modes^ of Model A (for £ ^ 2): since the 
lowest frequency mode has ujo2 — 0.29, it immediately 
follows that for 6 > 14 the w-mode frequencies can't be 
found in the Fourier spectrum. This argument is con- 
firmed by the analysis of the energy spectra, that are 
depicted in Fig. [51 The frequency distribution is consis- 



^ These numbers have been computed by a frequency domain 
code whose characteristics and performances are described in 

Refs. SEil,!!!,!!! 



tent with the value uj^^^^^ — 3\/2/b and thus the w-modes 
can be excited only for b < 14.4. Note that the differ- 
ent amplitudes of the spectra in Fig. [2] are due to the 
convention used for the normalization of the initial data. 

The same argument holds for the Black Hole. Since 
we have M = 1.4 the fundamental QNMs frequency is 
0.2669 (See Tab. HI] built from Table 1 of [ll); as a 
result, one must have b < 15.9 to trigger the fundamental 
mode (that dominates the signal) although the overtones 
(that have lower frequencies) can already be present in 
the waveform. In any case, w^ax smaller than the 
4th overtone only and, due to the correspondingly large 
damping time, this is not expected to give a recognizable 
signature in the waveform. 

On the basis of these considerations, we can summarize 
our results by saying that, for our M = 1.4M0 models, 
when b > 16, the incoming pulse is totally unaffected by 
the short-range structure of the object and the signals 
backscattered by any of the stars and by the Black Hole 
are identical in practice. This information, deduced by 
inspecting the waveforms, can be synthesized by compar- 
ing, as a function of b and for a fixed £, the energy released 
by the star (Sns) and by the Black Hole (Ebh) computed 
from Eq. (|2.7p . Figure [3] exhibits the ratio Ebh/E-^s for 
f = 2 and ^ = 3 (the latter for Model A only). This 
quantity decreases with b because (see Ref. [l[) for small 
b the Black Hole, contrarily to the star, partly absorbs 
and partly reflects the incoming radiation. On the other 
hand, the ratio tends to one for b > 16, in good numeri- 
cal agreement with the value of the threshold, needed to 
excite the quasi-normal modes, that we estimated above. 
Notice that the saturation to one for £ = 3 occurs for val- 
ues of b smaller than for £ = 2. This is expected: in fact, 
the QNMs frequencies increase with £ and thus one needs 
narrower b (and thus a larger wfjiax) to trigger space-time 
vibrations. 



TABLE II: The first four frequencies Un2 and damping times 
T„2 of ui-modes (for ^ = 2) of Model A: they have been 
computed by means of a frequency domain code described 
in Ref. HO, [23, • The third and fourth column of the table 
list the corresponding complex frequencies cj„2 — ian2 in our 
standard units. We have ujni = l-nVn^M^G j (? . 



n i/„2 [Hz] 


rn2 [ms] 


<^n2 


OLn2 


9497 


32.64 


0.29393 


0.15091 


1 16724 


20.65 


0.5176 


0.23853 


2 24277 


17.21 


0.75136 


0.28621 


3 32245 


15.43 


0.99796 


0.31923 


TABLE HI: The first four complex £ 
ujn2 — iQn2 of a M = 1.4 Black Hole 
(Derived from the values published in 


= 2 QNMs frequencies 
in our standard units 
Table 1 of [sj). 


n Un2 [Hz] 


T„2 [/is] 




an2 


8624 


77.52 


0.2669 


0.0635 


1 8002 


25.18 


0.2477 


0.1956 


2 6948 


14.42 


0.2150 


0.3416 


3 5804 


9.78 


0.1796 


0.5037 



B. Identification of the lu-modes 

We conclude this section by discussing the possibility 
of identifying unambiguously the presence of w-modes in 
the waveforms and in the corresponding energy spectrum. 
Ideally, one would like to find precise answers to the fol- 
lowing points: (i) understand which part of the waveform 
can be written as a superposition of w-modes; (ii) how 
many modes one should expect to be excited and (iii) 
how does this depend on h. 

Although these questions have been widely investi- 
gated in the past (see for example Chapter 4 of [1], 
Ref. [Q] and references therein), still they have not been 
exhaustively answered in the literature. The major con- 
ceptual problems underlying this difhculty are (i) the fact 
that the quasinormal-modes sets are not complete and 
(ii) the so called time shift problem. The former is in- 
trinsic in the definition of the quasinormal modes and 
prevents, in fact, to associate an energy to each excita- 
tion mode. The latter is related to the exponential decay 
of the quasinormal modes and it implies that, if the same 
signal occurs at a later time, the magnitudes of the modes 
will be larger with respect to that of the same signal oc- 
curred at an earlier time. As a consequence, the use of 
the magnitude of the amplitudes C„ (see Eq. 13.11 below) 
is not a good measure of the excitation of the quasinor- 
mal modes. We refer to the review of Nollert [35] for a 
thorough discussion of such problems. 

Beside these conceptual difficulties, from the practi- 
cal point of view it is however important to extract as 
much as information as possible about the quasi-normal 
modes by analyzing the ringing phase of the signal. Two 
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FIG. 4: Fits of the ring-down part of the waveform with the 
fundamental (n = 0) space-time mode for a Black Hole ex- 
cited by a 6 = 2 Gaussian pulse. We show the waveform 
^I'j"'' it) and its absolute value on a logarithmic scale to high- 
light the differences with the fit. 

complementary methods can be used to obtain such im- 
portant knowledge. On the one hand, one can imple- 
ment the Fourier analysis, namely looking at the energy 
Fourier spectrum in the frequenc y ra nge where w-modes 
are expected (see e.g. Refs. [H,!!!, [llj). On the other 
hand, one can perform a "fit analysis" . In this case, it is 
assumed that, on a given interval Au = [m;, Uf], the wave- 
form can be written as a superposition of n exponentially 
damped sinusoids, the quasi-normal modes expansion: 

= ^ = X! ^" cos(tj„u-f 0„) exp(-Q;„u) , (3.1) 

11=0 n=0 

of frequency w„ and damping time 1/q;„, that are, a pri- 
ori, unknown [we omit henceforth the index ^ since in the 
following we will be focusing only on the 1 = 2 modes]. 
Using a non-linear fit procedure one can estimate the 
values of (cJn, q;„, Cn 4>n) fi'om the waveform. We per- 
form this analysis by means of a modified least-square 
Prony method (see e.g. the discussion of Ref. 34]) to fit 
the waveforms. A feedback on the reliability of our fit 
procedure is done by comparing the values of frequency 
and damping time, LOni and Q;„f , obtained by the fit with 
those of Table [III and Table [Till that we assume to be the 
correct ones. 

The typical outcome of the fit analysis, using only the 
fundamental mode (n = 0), are shown in Fig. [4|for the 
Black Hole with 6 = 2 and in Fig. [Sj for the star Model A 
with 5 = 2 (top panel) and 6 = 8 (bottom panel) . When 
6 = 2, for which the largest space-time mode excitation 
is expected, for both the star and the Black Hole the fits 
show excellent agreement with the numerical waveform at 
early times, that progressively worsen due to the power- 
law tail contribution. The reliability of the procedure is 
confirmed by the values of ujq and that we obtain from 
the fit. 
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FIG. 5: Fits of the ring-down part of the waveform with the 
fundamental (n — 0) space-time mode for the stellar Model A 
exited hy b — 2 (upper panels) and a 6 = 8 (bottom panels) 
Gaussian pulse; For each value of 6 we show the waveform 
$2°'(i) a-nd its absolute value on a logarithmic scale to high- 
light the differences with the fit. 



For the Black Hole, we have uiq = 0.2660 and ao = 
0.0631, which differ of respectively 0.3% and 0.6% from 
the "exact" values of Table IIIII We can thus conclude 
that the fundamental mode is essentially the only mode 
excited for b — 2. For the star. Model A, we obtain 
ujo = 0.2739 and ao — 0.1636 and they differ respec- 
tively of 7% and 8% from the "exact" values of Table [Tll 
Since the damping time of the fundamental star w-mode 
is generically smaller than that of an equal mass Black 
Hole, the ringing is shorter and it is more difficult to 
obtain precise quantitative statements. In this case we 
tried to include more modes in the template (j3.ip used 
for the fit in order to precisely quantify the real contri- 
bution due to the presence of overtones in the signal. 
Unfortunately, in this case the fit procedure seems badly 
conditioned and we could not obtain a sensible feedback 
of the frequencies even if we clearly obtain (having more 
adjustable parameters) a better fit. 



BH b=2 




120 125 130 135 140 145 



u. 



Model A b=2 
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FIG. 6: The residual 7^ = 1-9 (See Eq. E^]) of the fits of the 
waveform of the response to a fe = 2 Gaussian pulse of a Black 
Hole (upper panel) and stellar Model A (bottom panel) as a 
function of the initial time (ui) of the fitting window around 
its best values that it is Ui = 129 for a Black Hole and Ui = 131 
for Model A. 



We have found that the choice of the time window to 
perform this analysis has a strong influence on the re- 
sult of the fit. This choice is delicate and it is related 
to the aforementioned problem of the time shift. Ideally, 
the window should start with the ring-down (i.e., at the 
end of the burst) and it must be both sufficiently narrow, 
in order not to be influenced by the non-oscillatory tail, 
and sufficiently extended to include all the relevant in- 
formation. There are no theoretical ways to predict or 
estimate the correct window, but some systematic proce- 
dures have actually been explored in the literature. We 
decided to use a method very similar to the one discussed 
in details in Ref. [13]: it consists in setting Uf at the end 
of the oscillatory phase, which is clearly identifiable in 
a logarithmic plot, and choosing the initial time of the 
window Ui such as to minimize the difference between the 
real data (vp^ata^ ^^^^ ^^le waveform synthesized from the 

results of the fitting procedure (^'^'). This difference is 
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BH b=2 



Model A b=2 




0.262 



BHb=2 





138 



FIG. 7: Determination of the best window for the fit of the Black Hole waveform and Model A (6 = 2). The initial time Mi is 
chosen so to minimize the residual TZ between the data and the fit. For the Black Hole we obtain Ui = 129 (u{ = 254), while 
for Model A we have Ui = 131 {u{ — 174). The horizontal line indicates the "exact" values of the considered model reported in 
Table HI] and Table Hill 



estimated by means of the following "scalar product" 



fit\ 



(3.2) 



whose result 8 is a value in the interval [0, 1] that it is 
exactly one when the two time series are identical (perfect 
fit). Fig.[S]shows such a determination for the Black Hole 
(top panel) and Model A (bottom panel) with b — 2: 
both curves exhibit a clear minimum of the quantity TZ = 
1 — 8 at, respectively, Mi = 129 and Ui — 131. The time 
window extends to U{ = 254 (Black Hole) and Uf = 174 
(Model A), respectively. 

As can be seen in Fig. [7] one has that even a small 
change of initial time Mi of the window used produces sen- 
sible variation of the estimated values of the frequency 
and the damping time of the fundamental mode. How- 
ever, it should be noticed that the estimated values ob- 
tained for the best window are those in best agreement 
with the expected values. 



We finally repeat the analysis for the datasets relative 
to wider Gaussian pulses. Focusing on the representa- 
tive 6 = 8 case, we find essentially the same picture with 
two main differences. First, the "global" quality of the 
fit, given by TZ is less good than in the 6 = 2 case; this 
is particularly evident for the star, where the fitted fre- 
quencies differ from the "exact" ones by more than 10%. 
Second, the fitting window becomes narrower and nar- 
rower as 6 is increased (see column seven of Table lIVp , 
thereby the the fit analysis quickly becomes meaningless. 



For the reasons outlined above we have found this pro- 
cedure not as effective as we hoped. We think that the 
problem of the unambiguous determination of the "right" 
time interval for the fit and of the presence and quantifi- 
cation the overtones in numerical data deserves further 
considerations. 
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TABLE IV: The results of the fit of the fundamental mode for 
the Black Hole and for the stellar Model A in the response 
to Gaussian pulse with b = 2 and 6 = 8 for the best fit 
windows [mi, itf] determined using the minimum of the residual 
TZ = 1 — Q criteria (see Fig. [6]). The reported error refer 
to relative difference between the fitted values and reference 
values reported in Table HTl and Table Hill 



Model b LOO Sujo % ao Sao % [ui,Ui] TZ 



A 


2 


0.2739 


7 


0.1636 


8 


[131,174] 


2 X 10"* 


A 


8 


0.3396 


15 


0.1302 


13 


[131,154] 


1 X 10~^ 


BH 


2 


0.2660 


0.3 


0.0631 


0.6 


[129,254] 


4 X 10"" 


BH 


8 


0.2614 


0.2 


0.0591 


0.7 


[129,229] 


5 X 10""^ 



IV. CONCLUSIONS 

In this paper we have studied numerically the scatter- 
ing of odd-parity Gaussian pulses of gravitational radia- 
tion off relativistic stars and Black Holes. We have found 
that the excitation of w-modes and black hole QNMs oc- 
curs basically in the same way for both objects: pulses of 
small b (high frequencies) can trigger the w-modes, while 
for large b (low frequencies) one can only find curvature 
backscattering effects and non-oscillatory tails. When w- 
modes are present, we have shown that both frequency- 
domain (energy spectrum) and time-domain (fit to a su- 
perposition of w-modes) analysis are useful to understand 
the mode content of the waveforms; however, our study 
also indicates that it is difficult to single out precisely the 
contribution of each mode, since the fundamental mode 
always dominates the signal and a clear identification of 
the overtones is lacking in the case of excitation induced 
by the scattering of odd-parity Gaussian pulses of gravi- 
tational radiation. 



The inspiring idea of this paper was to understand the 
origin of the dynamical excitation of w-modes in a sim- 
ple, but rather general, setting, where it is possible to 
do many, quick and controllable high-accuracy numeri- 
cal simulations with a tunable "source" . Our expecta- 
tion is that the main features of the process of lu-mode 
excitation (i.e., its dependence on the intrinsic frequency 
content of the initial data) that we have highlighted in 
the perturbative regime are sufficiently robust to survive 
qualitatively also in the full 3D-nonlinear case that needs 
the solution of the full set of Einstein equations. An en- 
couraging motivation for hoping so is that in other phys- 
ical settings, like for example the merger of two equal 
mass Black Holes, the perturbative analysis of the early 
70s Q in the extreme mass ratio limit was able to single 
out the important physical elements (i.e., the presence of 
QNMs), sketching a picture that has been later refined, 
but substantially confirmed, by numerical relativity sim- 
ulations ^36,, .31,' 38, 39, 40]. 
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